Programmable purification of type-I polarization-entanglement 
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We suggest and demonstrate a scheme to compensate spatial and spectral decoherence effects in 
the generation of polarization entangled states by type-I parametric downconversion. In our device 
a programmable spatial light modulator imposes a polarization dependent phase-shift on different 
spatial sections of the overall downconversion output and this effect is exploited to realize an effective 
purification technique for polarization entanglement. 
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Entanglement is a key resource for several quantum 
communication and quantum computation protocols [l|. 
Entangled states of polarization qubit have been real- 
ized exploiting spontaneous parametric down conversion 
(SPDC) in nonlinear crystals |2| where a pump pho- 
ton generates a couple of daughter photons called sig- 
nal and idler. Brightness and purity of SPDC-based 
entanglement sources are limited by decoherence. Ba- 
sically, there are two effects that degrade type-I polar- 
ization entanglement. The first is related to the source 
coherence time and the second to the emission angle- and 
frequency-dependent phase of the SPDC mechanism it- 
self. The effect of a limited coherence of the source may 
be reduced using a temporal precompensation technique 
[3,|j| whereas entanglement purity is usually recovered by 
strong filtering in the angular and spectral SPDC distri- 
bution (e.g. narrow irises and spectral filters). Recently 
very bright source has been realized compensating the 
angular-phase using suitably prepared birefringent crys- 
tals [5]. 

In this Letter we suggest and demonstrate a emission- 
angle and emission-frequency compensation technique 
based on the use of a programmable spatial light mod- 
ulator (ID SLM). This is a crystal liquid phase mask 
(64 x 10mm) divided in 640 horizontal pixels, each wide 
d = 100/im and with the liquid crystal lOfim deep. The 
SLM is set on the path of signal and idler at a distance 
D — 500mm from the two generating crystals. Driven 
by a voltage the liquid crystal orientation in correspon- 
dence of a certain pixel changes. The photons with an 
horizontal polarization feel an extraordinary index of re- 
fraction depending on the orientation, and this introduce 
a phase-shift between the horizontal and the vertical po- 
larizations. Since each pixel is driven independently it is 
possible to introduce a phase function dependent on the 
position on the SLM, i.e. on the SPDC generation an- 
gles. This kind of devices have been already used as am- 
plitude modulators Q, as well as diffractive elements to 
operate on orbital angular momentum hi and for phase- 
modulation at the single-photon level Q. The value of 
our technique is twofold. On the one hand, our filtering 
of downconverted photons is drastically reduced, not only 
in the angular distribution but also in frequency. On the 
other hand, programmability make it easily adjustable 
for many implementations. 



In our scheme (see Fig. [T]) the polarization entangled 
state is generated by two adjacent type-I BBO crystals 
0) H) 13 > oriented with their optical axes aligned in per- 
pendicular planes and pumped by a 405nm CW laser 
diode (Newport LQC405 - 40P). The part of the pump 
horizontally polarized generates couples of vertically po- 
larized photons in the first crystal, while the vertical 
pump polarization generates couples horizontally polar- 
ized in the second one. An half-wave plate (HWP) is 
used to set the pump polarization at 45° in order to bal- 
ance the probabilities to generate vertical or horizontal 
couples of photons. 
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FIG. 1: (Color online) Experimental setup. Polarization 
entangled states are generated by downconversion in type-I 
BBO crystals pumped @405nm by a Newport LQC405 - 40P 
CW laser diode. Then, a spatial section of the output cones 
passes through the SLM which provides purification of polar- 
ization entanglement. After SLM there are (on both paths) 
a slit, filters and a lens that make the image of the pump 
spot on the crystal into the core of a multimode optical fiber 
(MMF). For tomographic reconstruction we insert a quarter- 
wave plate, a half-wave plate and a polarizer and for the op- 
timization of the phase functions only the polarizers. 



Taking care of the coherent spatial overlap of the 
two emission cones the state outgoing the two crys- 
tals may be written as \$>(oj p ,uj s ,0)) = l/\/2{\HH) + 
exp{up(uj p , uj s , 0)}\VV)], where cu p and ui s are the shift 
from the pump and signal central frequencies w^ and 
Wj = ujp/2, and 9 is the shift from the central signal 
angle 0$. Upon expanding all the contributions to the 
optical paths to first order we obtain the phase function 
<£>(w p , uj s , 0) = ipa + aLuj p + /3Luj s — SL6, where L is the 
crystals length and 4>q is a constant term depending on 
the central frequencies and angles. The phase term aLuj p 



accounts for the delay time between horizontal and ver- 
tical downconverted photons. The last two terms appear 
because the photons generated in the first crystal must 
traverse the second one. These terms may be understood 
by considering the conservation of the transverse momen- 
tum. Taking the signal at the fixed angle 80 + 9, for differ- 
ent lu s the idler sweeps different angles 8' +8'. This means 
a different optical path and explain the phase term /3Llo s . 
Likewise fixing u B > a positive variation of 8 correspond to 
a negative variation of 8' and this introduces an optical 
path dependent on 8, i.e. the phase shift SL8. The con- 
servation of the transverse momentum implies that the 
signal and idler generation angles 9q + 8 and 8' + 8' are 
not independent. We write 8' = —d + -fuj s + j'ujp. Within 
the pump spectral width the dependence on uj p is negligi- 
ble and thus we have 8' ~ — 8 + -fuj s , with 7 = 88' /dco s = 
-0.045 x 10~ 15 rad/Hz = 1.294 x 10~ 4 rad/nm. Using 
this relation and 6 — 2/3/7, which can be proved analyti- 
cally, the phase function may be written as <p(w p , 8, 8') = 
0o + aLuj p - fihLQ + (3/jL8'. We exploit the SLM to 
add the phase functions 4>(8) and <j>'(9') to the signal and 
the idler respectively. This introduces a polarization de- 
pendent shift on the \HH) component of \^(uj p ,uj s ,8)) 
which, together with temporal delay compensation on 
the pump, allows one to achieve complete purification. 
The overall output state may be written as 
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where A8 is the angular acceptance due to the pres- 
ence of the slit. A(cjp) is the complex amplitude spec- 
trum of the pump laser and / = sinc(l/2AkuL), where 
Afc|| = Ak\\(uip,ui s (d,8'),9) is the spectral/angular am- 
plitude of the downconverted photons [9J . We stress that 
first order effects, i.e. frequency- and angle-dependent 
phase-shifts at the phase mask are negligible due to the 
small thickness (= 10 urn) of the mask. This definitely 
show the programmable nature of our device. 

In order to assess our model and verify the approxima- 
tion 9' ~ — 9 + juj s we have measured the angular dis- 
tribution of the horizontal component of the idler upon 
setting the signal at the fixed angle 9q + 9. The angu- 
lar resolution is A9 ~ 1.2 mr ad, corresponding to a slit 
aperture of 1 mm. We verify the absence of diffraction 
effects using a CCD camera and measure the angular 
distribution of coincidences and single counts using two 
set of filters set in front of each detector. We evaluate 
numerically the distribution of coincidences Pk(9o,9' ) ex 
J Ae d8j Ae d8'F?\f(0, d + 9,9' + 9')\ 2 and single counts 

p k (8 ,8> ) ex J Ae>>1 d8 j Ae d6>F i \f(o,e + e,e + e')\ 2 (i = 

1,2) as obtained omitting integration over pump fre- 
quencies (it can be shown that f(ui p ,8o + 9,9' Q + 9') m 
/(0, 9o+9, 8' +8') for crystals length L < lmm) and using 
Fi = Fbp, a 10 nm bandpass filters and F 2 = F£ P x QE, 
i.e. a filter built by two Fpp longpass filters (cut-on 
wavelenght 715 nm) and the quantum efficiency QE of 



each detector (taken from the specification of the photo- 
diode C309215). In Fig. [5] we compare the experimental 
results (points) with the theoretical distributions (solid 
lines): for panels (a,b,c,d) 8q = 3°, whereas for panel (e) 
the signal is set at 8q ± 8, with 8 — 1.8 mrad. In Fig. [2] 
(a) and (e) the experimental data are deconvolved with 
the function of the slit [10| • As it is apparent from the 
plot we have very good agreement with the theoretical 
expectations. The mismatch in the right part of panel 
(d) is due to the finite size of the phase mask. 
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FIG. 2: (Color online) Angular distribution of the idler 
(around 9' — —3°) with signal angle set at #o + 8. Blue 
points are the experimental data; the red solid line the the- 
oretical prediction, a) Coincidences per second with Fbp, 
9 = 0; b) Coincidences per second with 2 x Flp, 9 — 0; c) 
Single counts per second with Fbp, 9 = 0; d) Single counts 
per second with 2 x Flp, 9 — 0; e) Coincidence per second 
with Fbp, 9 = ±1.8mrad (blue circle/red square). 

For the purification procedure we set the slits at 
A8 ~ Q.bmrad and use, on both channel, the two filters 
F L p. This allows us to collect the downconverted pho- 
tons within a wide spectrum and angular distribution. To 
collect as many photons as possible we make the imaging 
of the pump spot on the crystals (~ 1.5 mm) into the op- 
tical fibers core (diameter of 62.5/xto) using the coupler 
lenses; the coincidence counts are about 100/s, basically 
due to low quantum efficiencies of our home made de- 
tectors (< 10%). Notice that the angular acceptance Ad 
acts as a lOOnm bandpass spectral filter for the down 
converted photons. The delay time between the photons 
may be compensated upon the introduction of a proper 
combination of birefringent crystals on the pump path, 
as already demonstrated in @, |j| (see Fig. []}. Upon set- 



ting <t>'{8') = P/jL9' + e' and (f)(8) = -/3/-yL8 + e, with 
e' + e = 4>q, it is possible to compensate all the remain- 
ing phase terms in ip(u p ,8,8') and to achieve complete 
purification of the state. With this choice we may gen- 
erate polarization entangled state of the form (JXJ) with- 
out any phase term. It is worth noting that the SLM 
also replaces the birefringent plate used for the opti- 



mal generation of photon pairs |2f. Using the SLM we 
set the phase functions 4>{x) — a,2(x — x C 2) + 62, and 
4>'(x) — a\{x — x c \) + 61, where x is the pixel number, 
x — x C 2 = -j 6* and a; — x c i = ~r9', X c i and £ C 2 are the 
central pixels on idler and signal, i.e. the pixels corre- 
sponding to the central angles 9' and 9$. The values of 
the parameters ai, 61, 0,2 and 62 has been optimized upon 
inserting two polarizers set at a.\ = 45° and a2 = —45° 
in front of the couplers and then searching for the min- 
ima in the coincidence counts, corresponding to the val- 
ues of 6i ; 2 compensating the constant phase difference 
4>o and ai_2 removing the angular dependence on 9 and 
& . For our configuration a\ = —02 = f3Ld/"fD ~ —0.05, 
b\ + b 2 = 0o- 
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FIG. 3: (Color online) Coincidence counts on a time window 
of 30s with the polarizers in front of the couplers set at 45° 
and —45°. (a) Coincidences as a function of 61,2 (blue/red) 
with &2,i = and optimal ai,2; (b) Coincidences as a function 
of ai — —a,2 with optimal 61 and 62 = 0. 

In Fig. [3£a) we report the coincidence counts on a time 
window equal to 30s as a function of 61 (62) (blue/red) 
with 62 = (b\ — 0) and with a\ t 2 set to their optimal 
values. In Fig. [3Jb) we report the coincidence counts on a 
time window of 30s as a function of a\ = —a 2 with 62 = 
and bi = ())q. The agreement with the theoretical model 
is excellent. In turn, the purification of the state works as 
follows: starting from a visibility equal to 0.423 ± 0.016 
we achieve 0.616 ± 0.012 after the delay compensation 
with the crystals and 0.886± 0.012 after the spatial com- 
pensation with the SLM. The residual lack of visibility 
is due to low spatial coherence of the pump, which is 
spatially multimode. Upon adding a constant phase £ to 
the term b\ or 62 we may also generate a wide zoology of 
entangled states. With £ = 0, it we have, in ideal condi- 
tions, the Bell states l^*) whereas with ( — ir/2 we may 
generate the state |$ w/2 ) = l/y/2[\HH) + i\VV)\. 

In order to characterize the output state and to check 
the effects of the decoherence processes, we have per- 
formed state reconstruction by polarization qubit to- 
mography. The procedure goes as follows: we mea- 
sure a suitable set of independent two-qubit projectors 



[111 [l2J corresponding to different combinations of polar- 
izers and phase-shifters and then the density matrix is 
reconstructed from the experimental probabilities using 
maximum-likelihood reconstruction of two-qubit states. 
The tomographic measurements are obtained by insert- 
ing a quarter-wave plate, a half-wave plate and a polar- 
izer. The results are summarized in Fig. 2] Fidelities 
of the purified polarization states |$) = |$ + ), |$~) and 
l^-ff/2) are F ~ 0.90 ± 0.01, where uncertainty has been 
evaluated by propagating the errors on the determination 
of matrix elements by tomography. In order to achieve 
this precision we have employed a long acquisition time 
(~ 60s) thus also demonstrating the overall stability of 
our scheme. The limited value of the fidelities is basically 
due to spatial incoherence of the pump and the incom- 
plete temporal delay compensation. 




FIG. 4: (Color online) Characterization of the output state. 
We report the tomographic reconstruction of the purified 
polarization entangled state (a)|$ + )($ + | (real part), (b) 
|$~)($~| (real part), (c) real and (d) imaginary part of 

|$7r/2)<$,r/2|. 



In conclusions, we have suggested and implemented a 
scheme to compensate spatial and spectral decoherence 
effects of type-I SPDC-based entangled states. In our 
device a programmable spatial light modulator acts on 
different spatial sections of the overall downconversion 
output and provides polarization entanglement purifica- 
tion. Our device allows the effective generation and pu- 
rification of type-I polarization entanglement and paves 
way for quantum information processing with low-cost, 
low-coherence, low-power source of entanglement |13| . 



[1] I. L. Chuang and M. A. Nielsen, Quantum Information 
and Quantum Computation, Cambridge University Press 
(Cambridge, UK, 2000). 

[2] P.G. Kwiat, E. Waks, A. G. White, I. Appelbaum, and 
PH. Eberhard, Phys. Rev. A 60, R773 (1999). 

[3] Y. Nambu, K. Usami, Y. Tsuda, K. Matsumoto, and K. 



Nakamura, Phys. Rev. A 66, 033816 (2002). 
[4] S. Cialdi, F. Castelli, I. Boscolo, M.G.A. Paris, Appl. 

Opt. 47, 1832 (2008). 
[5] R. Rangarajan, M. Goggin, and P. Kwiat, Opt. Expr. 17, 

18920 (2009). 
[6] G. Lima, A. Vargas, L. Neves, R. Guzman, and C. Saave- 



dra, Opt. Expr. 17, 10688 (2009). [11] K. Banaszek, G.M. DAriano, M.G.A. Paris, and M.F. 
[7] E. Yao, S. Franke-Arnold, J. Courtial, M.J. Padgett, and Sacchi, Phys. Rev. A 61 010304(R) (1999). 

S.M. Barnett, Opt. Expr. 14, 13089 (2006). [12] D.F.V. James, PG. Kwiat, W.J. Munro, and A.G. 
[8] S. Cialdi, I. Boscolo, F. Castelli and V. Petrillo, New J. White, Phys. Rev. A 64, 052312 (2001). 

Phys. 11, 023036 (2009). [13] S. Cialdi, D. Brivio, M.G.A. Paris, Phys. Rev. A 81, 
[9] S. Cialdi, F. Castelli and M.G.A. Paris, J. Mod. Opt. 56, 042322 (2010). 

215 (2009). 
[10] S. M. Riad, Proc. IEEE 74 , 82 (1986). 



